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A model-independent search for deviations from the Standard Model prediction is performed 
in e~p collisions at HERA II using HI data recorded during the years 2005 and 2006, cor- 
responding to an integrated luminosity of 159 pb~^. All event topologies involving isolated 
electrons, photons, muons, neutrinos and jets with high transverse momenta are investigated 
in a single analysis. Events are assigned to exclusive classes according to their final state. A 
statistical algorithm is used to search for deviations from the Standard Model in the distribu- 
tions of the scalar sum of transverse momenta or invariant mass of final state particles and to 
quantify their significance. A good agreement with the Standard Model prediction is observed 
in most of the event classes. No significant deviation is found in the phase-space and event 
topologies covered by this analysis. 



1 Introduction 

At HERA electrons and protons collide at a centre-of-mass energy of up to 319 GeV. These 
higti-energy electron-proton interactions provide a testing ground for the Standard Model (SM) 
complementary to e+e~ and pp scattering. 

The approach described in this papei^ closely follows the strategy of the previously published HI 
analysis using HERA I data . It consists of a comprehensive and generic search for deviations 
from the SM prediction at large transverse momenta. The analysis covers phase-space regions 
where the SM prediction is sufficiently precise to detect anomalies and does not rely on assump- 
tions concerning the characteristics of any SM extension. Using the complete HERA II e~p 
data sample, this is the first general search performed on a large data set from electron-proton 
collisions. 

° In this paper "electrons" refers to both electrons and positrons, if not otherwise stated. 
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Figure 1: The data and the SM expectation for all event classes with observed data events or a SM expectation 
greater than 1 event. The analysed data sample corresponds to an integrated luminosity of 159 pb~^. The error 
bands on the predictions include model uncertainties and experimental systematic errors added in quadrature. 



2 Data analysis 



The event sample studied consists of the full 2005-2006 HERA II e~p data set, correspond- 
ing to an integrated luminosity of 159 pb~^. All final states with at least two objects with 
Pt > 20 GeV in the polar angle range 10° < < 140° are investigated. Considered objects are 
electrons (e), photons (7), muons (/i), jets (j) and neutrinos (u) (or non-interacting particles). 
The identification criteria for each type of object are similar to those applied in the published 
HERA I analysis, ensuring an unambiguous identification while keeping high efficiencies. All 
objects are required to be isolated from each other by a minimum distance i? of 1 unit in the 
1] — (p plane. The events are classified into exclusive event classes according to the number and 
types of objects. This exclusive classification ensures a clear separation of the final states and 
allows an unambiguous statistical interpretation. 

As this analysis investigates all final state topologies of ep interactions at high transverse momen- 
tum, a precise and reliable estimate of all relevant HERA processes is needed. Hence, several 
Monte Carlo generators are used to generate a large number of events in all event classes, 
carefully avoiding double-counting of processes. The simulation contains the order as matrix 
elements for QCD processes, while second order a matrix elements are used to calculate QED 
processes. Additional jets are modelled using leading logarithmic parton showers as representa- 
tion of higher order QCD radiation. All processes are generated with a luminosity significantly 
higher than that of the data. 

The results of the data analysis are summarised in figure [U which presents the event yields 
subdivided into event classes for the data and SM expectation. All event classes with observed 
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Figure 2: The — log P values for the data event classes and the expected distribution from MC experiments as 
derived by investigating the ^ Pt distributions (left) and M^y[ distributions (right) with the search algorithm. 

data events or a SM prediction greater than 1 event are shown |. In each class, a good agreement 
between the number of observed data events and the SM prediction is seen. 
No data events are observed in the event classes i^-j-v and e-j-v where the largest discrepancy 
between the data and the SM prediction was found in the analysis of the HERA I data (see 
section |4]). Those classes correspond mainly to high Pt W production with subsequent leptonic 
decay, where deviations in the e~^p data continue to be observed^^l. The total SM expectation 
amounts to 1.2 it 0.2 and 2.5 it 0.8 in the ^-j-v and e-j-v classes, respectively. 

3 Search for deviations 

In order to quantify the level of agreement between the data and SM expectation and to identify 
regions of possible deviations, the invariant mass M^u and sum of transverse momenta Pt 
distributions of all reliable event classes are systematically investigated using the same search 
algorithm as developed for the previous publication. A region is defined as a sample of connected 
histogram bins, which have at least the size of twice the resolution of the observable. A statistical 
estimator p is defined to determine the region of most interest by calculating the probability that 
the SM expectation fluctuates upwards or downwards to the data. This estimator is derived from 
the convolution of a Poisson probability density function (pdf) to account for statistical errors 
with a Gaussian pdf to include systematic uncertainties. A possible sign of new physics is found 
if the expectation significantly disagrees with the data, and thus the region of most interest 
(greatest deviation) is given by the region having the smallest value, Pmin- This method finds 
narrow resonances and single outstanding events, as well as signals spread over large regions of 
phase-space in distributions of any shape. 

The fact that somewhere in the distribution a fluctuation with a value Pmin might occur is taken 
into account by calculating the probability P to observe a deviation with a p-value Pmin a-t any 
position in the distribution. Thus P is the central measure of significance of the deviation found. 
To determine P, hypothetical data histograms are diced according to the pdf of the expectation. 
The value of P is then defined as the fraction of hypothetical data histograms with a p^^^-valne 
smaller than the pj^in-value obtained from the data, and consequently the event class of most 
interest for a search is the one with the smallest P-value. 

The overall level of agreement between the data and SM expectation can be quantified further 
by taking into account the large number of event classes studied in this analysis. Among all 



The fjL — V event class is discarded from the present analysis. It is dominated by events in which a poorly 
reconstructed muon gives rise to missing transverse momentum, faking the neutrino signature. 



classes, there is some chance that small P values occur. This probability can be calculated by 
replacing all data distributions by hypothetical Monte Carlo (MC) distributions based on the SM 
expectation. The complete statistical algorithm is applied on this MC experiment, representing 
a single HERA experiment with an integrated luminosity of 159 pb~^. The expectation for the 
P values of the data is then given by the distribution of Psm values derived from many MC 
experiments. In the case that deviations arise from statistical or systematical fluctuations only, 
the distribution of P-values obtained from data and MC experiments are compatible. 

The results of the search for deviations between data and SM expectation are summarised in 
figure El Shown is the negative base-10 logarithm of the P values obtained from the real data 
compared to the expectation derived from a large set of MC experiments. The comparison is 
presented separately for the scans of the M^ii and J2 Pt distributions. All P values range from 
0.01 to 0.99, corresponding to event classes where no significant discrepancy between data and 
SM expectation is observed. These results are in good agreement with the expectation from MC 
experiments. 

Although data events are observed in the j-j-j-j and j-j-j-j-v event classes, no reliable P values 
can be calculated for these classes due to uncertainties of the SM prediction at highest M^n and 
^ Pt values^. These event classes are not considered to search for deviations from the SM in 
this extreme kinematic domain. Consequently, these event classes are not taken into account to 
determine the overall degree of agreement between the data and the SM. 

4 Comparison vi^ith HERA I results 

While good agreement in all event classes is observed between the data and the SM prediction 
in the present HERA II analysis, some discrepancy is found in the previously published general 
search on the HERA I data. Complementary to the pure e~p event sample studied here, the 
HERA I data set is largely dominated by positron-proton collisions. There the most significant 
deviation is found in the n-j-iy event class with P values of 0.01 and 0.001 for the scan of the 
Mg^ll and J2 Pt distributions, respectively. The global probability to find at least one event class 
with a P value smaller than that of the ^-j-v class in the HERA I data amounts to 28% for the 
Mgx\ and 3% for the Pt distributions. 

5 Conclusions 

The data collected with the HI experiment during the years 2005-2006 (HERA II) have been 
investigated for deviations from the SM prediction at high transverse momentum. All event 
topologies involving isolated electrons, photons, muons, neutrinos and jets are investigated in a 
single analysis. This is the first general search performed on a large set of data from electron- 
proton collisions. A good agreement between data and SM expectation is found in most event 
classes. In each event class the invariant mass and sum of transverse momenta distributions 
of particles have been systematically searched for deviations using a statistical algorithm. No 
significant deviation from the SM is observed in the phase-space and event topologies covered 
by this analysis. 
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